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Abstract 

We represent an algorithm ahowing one to construct new classes of partially inte- 
grable multidimensional nonlinear partial differential equations (PDEs) starting with the 
special type of solutions to the (l+l)-dimensional hierarchy of nonlinear PDEs lineariz- 
able by the matrix Hopf-Cole substitution (the Biirgers hierarchy). We derive examples of 
four-dimensional nonlinear matrix PDEs together with they scalar and three-dimensional 
reductions. Variants of the Kadomtsev-Petviashvili type and Korteweg-de Vries type 
equations are represented among them. Our algorithm is based on the combination of 
two Frobenius type reductions and special differential reduction imposed on the matrix 
fields of integrable PDEs. It is shown that the derived four-dimensional nonlinear PDEs 
admit arbitrary functions of two variables in their solution spaces which clarifies the in- 
tegrability degree of these PDEs. 

1 Introduction 

The problem of construction of such multidimensional nonlinear partial differential equations 
(PDEs) which are either completely integrable or, at least, possess a big manifold of particular 
solutions is very attractive problem of integrability theory. First of all, one should emphasize 
several remarkable classical works regarding the completely integrable models such as [T| (where 
the Korteweg-de Vries equation (KdV) has been first time studied), [21 [3] (where the dressing 
method for a big class of soliton and instanton equations has been formulated), [4j (where 
the so-called Sato approach to the integrability is applied to Kadomtsev-Petviashvili equation 
(KP)). A big class of the first order systems of quasilinear PDEs is integrated in [3, E] using 
the generalized hodograph method. However, most of the above integrable nonlinear PDEs are 
(2-1-1)- and/or (l-l-l)-dimensional, excepting the Selfdual-type PDEs [3 [8] and the equations 
assotiated with commuting vector fields [H [TOl [TTl [121 [ISl [E] • 

In this paper we use modification of the algorithm represented in [151 IB] allowing one 
to construct new multidimensional partially integrable nonlinear PDEs. It is shown in [16] 
that Frobenious reduction of the matrix fields of the nonlinear PDEs integrable either by the 
Hoph-Cole substitution [H] (C-integrable PDF [HI HH [201 [2ll [22l [23]) or by the method of 
characteristics ^24^125] (C^-integrable PDF) leads to one of two big classes of the nonlinear PDEs 
integrable by the inverse spectral transform method (ISTM) [2BII271 (S'-integrable PDEs 

[TS]). These classes are (a) soliton equations, such as KdV [Il[2n], the Nonlinear Schrodinger 
equation (NLS) [31], the Kadomtsev-Petviashvili equation (KP) [32], the Deavi-Stewartson 
equation (DS) [33], and (6) instanton equations, such as the Self-dual Yang-Mills equation 
(SDYM) [11 [8] . 
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The natural question is whether the Frobenius reduction (or its modification) can be used 
for construction new types of integrable (or at least partially integrable) systems starting with 
any known integrable system or this method works only for derivation of soliton and instanton 
PDEs from C- and C/i-integrable ones? 

At first glance the answer is negative. In fact, one can verify that Frobenius reduction 
applied to such matrix S'-integrable PDE as the GL(A^) SDYM, the A^-wave equation and the 
KP produces the same S'-integrable PDE. However, there is a method to generate new higher 
dimensional partially integrable systems of nonlinear PDEs using Frobenius type reduction after 
the appropriate differential reduction imposed on the matrix fields of the above S'-integrable 
nonlinear PDEs. 

Such combination of reductions has been already used in JM]. It is shown there that 
GL(A^) SDYM supplemented by the pair of reductions, namely, differential reduction relating 
certain blocks of the matrix field and Frobenius reduction of these blocks, produces a new five- 
dimensional system of matrix nonlinear PDEs (with three-dimensional solution space), whose 
scalar reduction results in the nonlinear PDE assotiated with commuting vector fields [131 E] • 

Following the strategy of ref.[31], we consider an algorithm for construction new partially 
integrable PDEs starting with the matrix KP (although this algorithm may be applied to any S'- 
integrable model). We will derive two representative of matrix systems whose scalar reductions 
yield the following four-dimensional equations: 

VyV:,yy = 0, (1) 

{vxd-\xx)}j = 0. (2) 
Another reduction, 

vt, = 0, (3) 
reduces these PDEs into the following three-dimensional ones 

1 3 

Vt3 - -^xxx - -^v^dy^v^x = 0. (5) 

Here we take r and as evolutionary parameters (times) while x, y and ^2 are taken as space 
parameters. Since the linear parts of eqs.([2]) and ([5]) coincide with the linear parts of the KP 
and the KdV respectively, eqs. ([2]) and may be treated as new variants of KP- and KdV-type 
equations respectively. The feature of these equations is that the derivative with respect to y 
appears only in the nonlinear parts. For this reason, these equations may be referred to as 
dispersionless ones. 

Remember, that the KP originates from the (l+l)-dimensional C-integrable Biirgers hier- 
archy due to the Frobenius reduction |TB] . Thus the complete set of transformations leading to 
eqs.([T]) and ([2]) is following (see also Fig.l). We start with the C-integrable Biirgers hierarchy 
of nonlinear PDEs with independent variables x and tn, n = 2,3, ... . Frobenius reduction of 
this hierarchy [16] yields the proper hierarchy of discrete chains of nonlinear PDEs, which is 
equivalent to the chains obtained in the Sato approach to the integrability of (2+l)-dimensional 
KP [1]. The later may be derived as an intermediate result of our algorithm after eliminating 



_1 _1 \ _1 

:t 2^2/*2 2^^0^xx j „^xxx ~l~ "^yyy^x 

' y 

3 1 3 / 

^xt'i T^t2t2 "T^xxxx ~l~ 2 yyyt2^x "^12^ xy 
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all extra fields using combination of the first and the second representatives of the constructed 
discrete hierarchy. Next, apply the differential reduction introducing one more independent 
variable y. This step yields a three-dimensional system of nonlinear PDEs, i.e. dimensionality 
coincides with that of KP. Finally, the Frobenius type reduction applied to the matrix fields of 
the latter nonlinear PDEs results in a four-dimensional system of matrix PDEs whose scalar 
versions yield eqs. (11121) . It will be shown that the derived four-dimensional nonlinear PDEs 
may not be completely integrated by our method because the available solution spaces to them 
are restricted. 
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Fig.l The chain of transformations from the (l+l)-dimensional C-integrable Burgers hierarchy to 
the four-dimensional PDEs ([T|) and Here N = 2NonoM, TT^*) are IuqM x 2noM matrix fields, 
w^^\ p^*) and 5^*) are ngM x ngM blocks, l2noM and Im are 2noM- and M-dimensional identity 
matrices respectively. Blocks p^^^ and q^^^ are defined by eqs. (j26p . Here A'^o, no and M are arbitrary 
positive integer parameters. 



This paper is organized as follows. In SecEJwe briefly recall some results of [I6] and derive 
the discrete chains of nonlinear PDEs produced by the (l-l-l)-dimensional Biirgers hierarchy 
(with independent variables x and t„, n = 2, 3) supplemented by the Frobenius reduction. Using 
a few equations of theses chains we eliminate all extra fields and derive the matrix KP. In Secl3] 
we suggest the special type differential reduction imposed on the blocks of the 2noM x 2noM 
matrix fields of the above chains. This reduction introduces a new independent variables y 
and r and allows one to generate ngM x uqM three-dimensional matrix PDEs. Frobenius 
type reduction of the above PDEs results in the four-dimensional systems of M x M matrix 
PDEs, SecJH Scalar versions of these PDEs result in eqs.([T]) and ([2]) which, in particular, may 
be reduced to eqs.Q and ([5]) respectively. Solution spaces to the nonlinear PDEs derived in 
Secsj3] and m will be studied in SecJSl The obstacles to the complete integrability of eqs.([I]), 
([2]), (jl]) and ([5]) as well as the obstacles to their integrability by the ISTM are briefly discussed 
in Secjni Conclusions are given in SecJTl 
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2 Relation between the (l+l)-dimensional C-integrable 
Burgers hierarchy and the matrix KP 



2.1 The (l+l)-dimensional C-integrable Burgers hierarchy 

Here we use the Burgers hierarchy as a simplest example of C-integrable hierarchies lineariz- 
able by the matrix Hopf-Cole substitution. Namely, let W be solution of the following linear 
algebraic matrix equation: 

X. = Wx, (6) 
Dt„X = d:x, A„x = Xt„+xA^"\ n = 2,3,.... (7) 

Here x ^^(^ ^ are 2NQnoM x2NQnoM matrix functions, A*^"^ are 2NQnoM x2NQnQM commuting 
constant matrices. To anticipate, the matrices A^") are introduced in order to establish the 
reductions eliminating derivatives with respect to tn from the nonlinear PDEs, see, for example, 
reduction ([3]). The compatibility conditions of eqs. ([6]) and (JTj) yield the linearizable Biirgers 
hierarchy. We will need only the first and the second representatives of this hierarchy below, 
i.e. n = 2, 3 in eqs. ([7]): 

Wt, - W,, - 2WJV = 0, (8) 

Wt.-w^^^-m^^w -?,w^{w, + w^) = Q. (9) 



2.2 The Frobenius reduction and assotiated chains of nonUnear PDEs 

Introduce the Frobenius reduction |16j : 



W 



02Mno hMno 
\ 02Afno 02Mno 



02Mno 02Mno 
02Mno 02Mno 

hMno 02Mno / 



(10) 



where Ij and Oj are J x J identity and zero matrices respectively, VF*^*-* are 2noM x 2noM 
matrix functions. Substituting eq. ffTUl) into eq.([n]) we obtain the following block structure of x: 



( 



X 



X 



(1) 



X 



{No) 



\ 



(11) 



where x'-*-' are 2noM x 2noM matrix functions. Thus, eqs. ([6]) and (j7j) may be written as follows: 



x^r\s) 



i=l 



= xt^ + x^"^A^'^^'"\ A^"^ = diag(A("^^), . . . , A^"^^")), 
m = 1, . . . , A^o, n = 2,3, . . . , 



(12) 

(13) 
(14) 
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where A*^"'™') are some 2noM x 2noM commuting constant matrices. The meaning of the integer 
parameters no and M is clarified in Fig.l. The compatibihty conditions of eqs. (fT2|) and (fT3|) 
(with n = 2, 3) yield the following discrete chains: 

W}^^ - Wi2^ - 2M^i"+') - 2M^i^)M^(") = 0, (15) 
Wt^s'^ - Wi2l - 3 (vrii^ + + Wi'^ (iy(i)iy(") + vr("+^) + Vl^j")) + (16) 

p|H2)p^(n) ^ j = Q, n = 1, . . . , A^q, 

lyW = 0, / > Nq. 

Alternatively, these chains may be derived substituting eq. ffTOj) into eqs.([8]) and ([9]). 
2.3 Matrix KP and its scalar reduction 

The matrix KP is represented by the system of three equations involving fields W^^\ i = 1,2, 3: 
eq. (fT5|) with n = 1,2 and eq. (fT6|) with n = 1. Eliminating P^(^) and W^'^'> from this system one 
gets the following nonlinear PDE for the field VF*^^-*: 

where square parenthesis mean matrix commutator. In the scalar case this equation reduces to 
the following one, u = W^^^: 

/ 1 3 2\ 3 

[ut;--u^xx--u^)^-2'^t2t2 = (18) 
which is the scalar potential KP. 



3 Differential reduction of the matrix KP 

Let matrices x^'''' have the following block structure: 



x^'^= ..(2j-i) ], j = i,...,iVo, (19) 



where \Ef(™) are uqM x noM matrix functions. We introduce y-dependence of vE^^"^) by the 
following second order PDE: 

£(0) _ ^(™) = a^i™) + i^^('") + /x^(™), m = l,...,2A^o, (20) 

Here a, v and ^ are n^M x n^M diagonal constant matrices. The block structure of 'y^^^ ( fT9l) 
suggests us the relevant block structure of W^^^: 

^''=(p(') go!)' ^■ = l'---'^o, (21) 
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where w^^\ v^^\ q^^^ and p^^^ are uqM x rioM matrix functions. Now matrix 2noM x 2noM 
equations (fT2!) may be written as two tt-qM x 2NQnoM equations: 



where 



£(1) := = ^ + , (22) 

i=l 

£(2) := ^^.^ = ^ + , (23) 



^ = (^W,...,^(27V°)). (24) 
The compatibihty condition of eqs. fl22|) and fl23|) . 



^f=^^'\ (25) 



3/ 



yields the expressions for p^^'' and g^-'-' in terms of v^^^ and tf^-^-': 

= ^(i) + ^(i)^ + ^a+i)^ + ^;(i)a«;(^'\ = v'y^^ + v^^v + + v^^\v^^ . (26) 

Thus only two blocks of W'^^^ are independent, i.e. w'^^^ and f '•■'^ Matrix equation (122|) may be 
considered as the uniquely solvable system of INq linear n^M x n^M matrix algebraic equations 
for the matrix functions v^'^^ and ly^*-', z = 1, . . . , A^o^ while eq. (!23|) is the consequence of eq. (l22!l . 
Eqs.([I3]) yield: 

Z),^^M = S^^^"^), A,,^^""^ = ^Ir^ + ^("^)a("''"), m = 1, . . . , 2Aro, (27) 
A("'^) = diag(A("'2^-i),A("'2^')), j = l,...,iVo, n = 2,3,..., 

where A'^"''"^ are n^M x noM commuting constant matrices. Eq. fl20|) allows us to introduce one 
more set of parameters r„ as follows: 

^M^^n^M^ m = l,...,2iVo, n = l,2,.... (28) 
Before proceeding further we select the following three equations out of system (I27|28p : 



Z),^^M = (29) 
^M^^M=^M^ (30) 
^^^^M ^ vpM. (31) 

We also assume that a, v and ;U are scalars, i.e. 

= -^noAf, V = T^O^noM-, ^J' = ^^O^noM■ (32) 

In order to derive the nonlinear PDEs for and we must consider the conditions providing 
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the compatibility of eq.( l22i) and eqs.( l29ll3Ti) : 

Dt,S'^'^ - Si'J = ^ (33) 

No No 

^ E['^^d;^^'Dy^ + i^r^r^'^ = 0, n = 1, 2, . . . , 

i=l 1=1 

f«-fS) = (34) 
^ E['^^d-^^'% + J2 Ei''^d-^^'^ = 0, n = 1, 2, . . . , 

i=l i=l 

A3^«-4i. = ^ (35) 
^ + ^ = 0, n = 1, 2, . . . . 

i=l i=l 

These equations generate the following chains of nonlinear PDEs for v^^^ and w^'^\ i = 1, . . . , Nq: 

e[''^ = t;«-^«-2i;(*+i)-2(i.ot;i') + «;«+i;«^;«)i;«-2^;^^^ (36) 
) = u;g) - «;« - 2t/7i^+i) - 2(4^) + - 2t;i')(«;i*) + ^^'""'^ + ^o^^^^) = 0(37) 

^ ^« _ ^« _ ^(m) _ ^^^W _ ^« _ (^^^(1) + ^« (38) 
i?f ) ^ «;« - w^l - ^) - - ^o^« - (..« + v^^^ + <).«)^.» - (39) 



i^f ^ - -ff - v^L - Se^^ - 3.r^) - 3 (wP +w^ + + + (.i^))^ + (40) 

^(i)(^(i))3 + ^^(^(1) + ^(1)^(1) + ^(1)^(1) + ^(2) + ^(1)^(1) + ^,^1)^(1) + 2vi'\v^'^f) + 

) ^ ^g) _ _ 3^(;+i) _ 3^(^+2) _ 3 (^^a) + ^(2) + ^(2)^(1) + ^(i)^(i) + ^a)^(i) + (41) 
+ ^(1)^(1) + ^(1)^(1) + (^(i))3 + ^(2) + ^(1) + ^(D) + ^(i)(^(i))2 + 

^,0(^(1)^(^+1) ^ ^(i)^W) ^ ^^(^(1)^(0 + ^(1)^(0 + ^» + (^;(i))2^;(0 + ^(^+1) + 
^,^^(1)^(0)^ _ 3t;(2)(t;(^+i) + + /zot;«) - 3«;(.i) + + ^;(i)(^;(^+i) + 



In addition, we must to take into account the compatibility condition of eqs.( l20l) and ( l22l) . 

1=0 

which gives us the following system of non-evolutionary nonlinear chains: 

E? := + 2«;« - + z/o4') + 2t;«t;« = 0, (43) 
■■= w^^-uow^;^ -w^J^ +2v^y'+^^ + 2fiov'^;^ +2v^^^w^'^ = 0, t = l,...,No. (44) 

One can show that the derived chains fl36ti40|H3P^ generate three compatible three-dimensional 
systems of nonlinear PDEs as follows: 

= 0, ^ =0, z = 0, 1, fields v^^\ w^^\ j = 1, 2, (45) 
independent variables t2,x,y, 
Ef^^ = 0, ^ =0, z = 0, 1, fields v^^\ w^^\ j = 1, 2, (46) 

independent variables T,x,y, 
Ef^^ = 0, ^ =0, i = 0,l, k = l,2, fields v^^\ w^^\ j = 1, 2, 3, (47) 

independent variables t3,x,y. 

We do not represent these equations explicitely since they are intermediate equations in our 
algorithm. 



4 The Frobenius type reduction, new four-dimensional 
systems of matrix nonlinear PDEs and their scalar re- 
ductions 

Three systems f l45M7j) represent three commuting flows with times t2, r and respectively. 
In this section we follow the strategy of ref. [TB] and show that hierarchy of nonlinear chains 
(1361141 supplemented by the Frobenius type reduction of the matrix fields f and ly*^*^ 

generates the hierarchy of four- dimensional systems of nonlinear PDEs, see eqs. (I65H68I) and text 
thereafter. 

This is possible due to the remarkable property of chains of nonlinear PDEs (1361141 |H3P^ . 
Namely, these chains admit the following Frobenious type reduction: 

yii) = [^m)^ fc,/ = l,...,no}, = A;, ^ = 1, . . . , no}, (48) 

v^''^^^ = 6kiv^''^^ + 5k{i+m(i))lM. 

= 4iw(*'') + 4(«+n2(i))/A/, i = 1, . . . , A^o, 

^(i;0 = = Oa/ , / > no, V ^, 

where tiq is an arbitrary positive integer parameter, ni{i) and ^2(2) are arbitrary positive integer 
functions of positive integer argument, f and are M x M matrix fields. In particular, 
if ni{i) = 712(2) = 1, then this reduction becomes Frobenius one pTO], which is shown in Fig.l. 
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Eqs.f l48p require the following block-structures of the functions ^^^^ and of the constant 
matrices A*-*''); 

^(0 = |^a;nm)^ n,m = l,...,no}, (49) 
X{i;0 = {a(*;';™), m = l,...,no}, l = l,...,2No, i = 2,3. 

Here vl/^';"-™-) are M x M matrix functions and A*^*'''*") are M x M commuting constant matrices. 
In turn, eq. (12^ reduces to the following one: 

A^o no 



, 1=1,.. . ,2iVo, n,m = l,...,no, 



i=i j=i 



while eqs.fl29l l3T]|20|) yield 

p^^^a;nm) _ ^a^H^ (51) 

^a;"'") = ^('^™), (52) 

p^^(Z;nm) _ _|_ ^{«;nm)^(i;/;m) ^ Z = 2, 3, 

/ = !,.. . ,2A^o, n,m = 1, . . . ,no. 

Then the chains of nonlinear PDEs 0361141 get the following block structures: 

Ei™) = A:,/,= l,...,no} = 0, n = l,2,3, (55) 

E^ = {El^'^6,u k,l,= l,...,no} = 0, 
m = 0, 1, i = 1, . . . , Nq. 

In particular, nonlinear chains fl36|37l) and fl43|44l) with Frobenious type reduction fB5]) generate 
the following chains with two discrete indexes: 

^(ii;0 ^ ^g;0_^^/)_2^C*+i;0_2(z/ot;(^^i) + t/;(i^i) + ^^^^^ (56) 
^(H;0 ^ _ ^(^0 _ 2^(^+i;0 _ 2{w^^^) + v^^'''^ v'^'-''^)w'^'-''^ - 2v^^''> (w'^^'^ + v^'^''^'^ +(57) 

Ef) := ^£0 + 2ii;f)-t;(^^nz^z;f) + 24^^i)^(^^') + 2t;f+"i(^ (58) 



■= ^if - ^o<''^ - w^^''^ + 2^;(^+i'') + 2fiovl'''^ + 2vl^'^^w^'''^ + 2v^''+''^'^^^ = 0, (59) 
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where i = 1, . . . , Nq, I = 1, . . . , tiq. Now our goal is to write the complete system of PDEs for 
some fields v^^'^^ and/or ty^"-'*") which would be independent on both parameters A^o and no (in 
the spirit of the Sato theory [1]). This may be done if, along with system (15611591) . one involves 
the discrete chains of PDEs generated by either eqs. (l38|39l) or eqs. (l40f41l) together with the 
Frobenious type reduction (HHI) . For instance, system (I38|39p yields 

j^m) ^ _ ^(,0 _ _ ^^^ir,l) _ _ ^^^^(m + ^(1;1) + ^(1;1) + (gg) 

^(l;l+ni(l))^(i;0 _^ ^(1;1) ^(l;/+ni (i)) _^ ^(l;«+ni (l)+ni (i)) _ ^{l;l+n2{i)) ^ 

^(l;l)(^(i+l;0 ^ ^(i;/) _^ _ ^(l;;+n2«) _ ^(l;«+n2W) _ (1)) ^(i;0 + 

^(l;l)^(l;«+n2(i)) _^ ^(l;/+ni (l)+n2 (i)) ) _ ^(l;/+ni (i+1)) _ = Q. 



The discrete chains generated by eqs. (140!) and pTl) are very cumbersome. However, it will 
be shown that we need only equation (1401) with i = 1 for derivation of the complete system 
of nonlinear PDEs. The later equation may be simplified after elimination v^"* and w^x'^ 
using eq. (l36i) with 2 = 1,2 and eq.(l371) with i = 1. This simplification is reasonable because, in 
any case, we are going to use chain (!56|) (produced by eqs. (136!) ) in order to derive the complete 
system of nonlinear PDEs. One gets in result: 

<^ - - 3^;^ - 6.Jf -Q(wf^ + + {v^^f + .ovf^y^ - (62) 



0. 



After the Frobenius reduction, this equation generates the following discrete chain: 



- - 3.S^ - 6.;f ) - 6(..if ) + .Jf ).(^^^) + (.(^^^))^ + .o.Jf - (63) 



6('u;(l;'+"l(l)) + ^(l;«+2ni(l)) ^ ^(l;l+ni(l))^(l;/) ^ ^(1;1) ^(l;/+„i (1)) ^ ^(l;l)^(l;«+ni(l)) ^ 



Having chains of equations fl56H61ll^^ we would like to exhibit the complete system of 
nonlinear PDEs independent on both Nq and no taking a few equations out of these chains. 
After some examination of these chains we conclude that it is enough to take the equations 
e[''^'^\ k = 1,2, e[^^'^^ and E^'^\ Introducing the new fields 

V = v'^'-''\ g = t;(l^l+"i«)+«;(l;l), (64) 

p = + + + ^(l;l+2ni(l)) ^ ^^^da+md)) ^ 
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we write these equations as follows: 

E^i^''^^ = Vt^ - v^^ - 2p^ - 2{uoV^ + qx + v^v)v - 2v^{q + Vy) = 0, (65) 

^'^^ = Vr - v^y -Py- VqV^ - qx- {i^oVy + Qy + Vx + Vyv)v - Vy{q + Vy) = 0, (66) 

-Sf ^'^^ := Avts - Vxxx - 3vxt2 - Qpt2 - 6 {qt2 +^'43^ + {^xf + v^Vt^v - (67) 
6f X [qx + v^y + vv^ + VqV^ - 6vt^ (q + Vy) - 6q^v^ = 0, 



e[^'^^ := Vyy + 2qy - Vx + jyoVy + '^VyV = 0. (68) 

Any two equations out of the system fl65ti67p supplemented by eq. (l68l) represent the complete 
system of matrix nonlinear PDEs for the fields v, p and q. Two examples of the scalar nonlinear 
PDEs are given below in Sec J4.1[ 

4.1 Examples of scalar nonlinear PDEs 

Example 1. Consider the system fl65|66|68l) . In the scalar case this system reduces to the 
following single scalar nonlinear PDE for the field v. 



{Er'%-^iEr''% +-iEr')xx = eq.©. (69) 

y 

Example 2. Consider the system fl65|67|68l) . In the scalar case this system reduces to the 
single equation for the field v as follows: 



^^r^ + ^5;^(^r^).)^-^(i?r^k eq.dSD. (70) 



5 Solutions to the derived nonlinear PDEs 
5.1 Solutions to nonlinear system (I45H47I) 



Solutions to systems of matrix nonlinear PDEs ( 14511471) or, more general, to discrete chains 
(1361141 |l43]|4^ are obtainable in terms of the functions \Ef("^) m = 1, . . . ,2A^o taken as solutions 
of linear system fl20ll29l l3T|) : 

= E / rf#gH9)e^^+'="^+(^'-^^'')*^-^('^'-^^"')*^+'^'''^^ (71) 
i=l ^ 

a,f3=l,...,noM, I = 1, . . . ,2No, 

where i^^liq) are the arbitrary scalar functions of single scalar variable, x = (x, ?y, ^2, ^3) is 
the list of all independent variables of the nonlinear PDEs, k^^\ i = 1,2, are roots of the 
characteristics equation assotiated with eq.( l20|) where a, y and /U are given by eqs.( l32|) : 

k"^ - q - uok - fio = 0, ^ (72) 
k^'^ = ^ (^0 + \/ vl + 4g + 4/io^ , A:(2) = \[vq- \J ^1 + 4.q + Aii^ . 
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Now, we can use eg. (1221) in order to find v'-*-' and w^^\ In particular, there is a big class of 
solutions in the form of rational functions of exponents, such as solitary wave solutions. We do 
not represent their explicite form since eqs. (I45ti47|l are intermediate equations in our method 
of solving eqs. ffT][2l) . 

Note nevertheless that the simplest nontrivial solution to eqs. (l45tH7|) corresponds to A^o = 2. 
Then eq. fl22p reduces to the following four matrix equations: 

4 

= + , m = 1, 2, 3, 4. (73) 

i=l 

These equations, in general, are uniquely solvable for the matrix fields t^^*^ and w^'^\ i = 1,2. 



The reductions to the lower dimensional PDEs. The described solution space admits 
the following reductions 

^« = u;« = 0, Vz, (74) 
where n is either 2 or 3. In fact, these reductions mean 

in eq. fITT]) . For instance, let n = 2. Then we must substitute the following expression into 
eq.dZI]): 

= - ^) + {l + ^) , (76) 

where ip^aft\ k = 1,2, are arbitrary constants. 

5.2 Solutions to system (I65H68I) and to scalar equations ([1]J2]) and 

dUSD 

In this case the fields v^^'^'^ and w^^'^'' are solutions to linear algebraic system (l50l) . More precisely, 
only part of system (!50l) is needed to define v^^'^^ and w^^'^K In fact, system (!50l) may be viewed 
as two subsystems. The first one corresponds to n = 1: 

TVo no 
i=l j=l 

1 = 1,.. .,2No , m = 1, . . . ,no. 

This is the system of 2iVono linear algebraic MxM matrix equations for the same number of the 
matrix fields v^^'''^ and w^^'''\ i = 1, . . . , Nq, j = 1, . . . , uq. Namely eqs. (1771) yield the functions 

^(l;l)^ ^(1;1)^ ^(2;1)^ ^(l;l+ni (1)) ^ ^(l;l+ni(l)) ^ (1)) ^ ^(l;l+2rM (1)) defining the SolutloUS V, p, q tO 

system fl65H68l) in accordance with formulae (I64p . The second subsystem corresponds to n > 1 
in dSU]): 

No 
i=0 

1 = 1,..., 2No, n,m = l,...,no, = 0, if, i < 0. 
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This subsystem expresses recursively the functions v|/(';"'^)^ n > 1, in terms of the functions 
^(«;im)_ rjj^g simplest case corresponds to nj{i) = 1, \fi,j (the Frobenius reduction). 
The functions v]/(';"'^) are solutions to system (I5m54p and may be written as follows: 

^iTH^) = E / dgC^^)(g)e^^'+'=^"^+(^^-4"''"^)*^+(''^-<'''"')*^-^^''^''- (79) 
a, (3 = 1, . . . ,M, 171 = 1,..., no, 1 = 1,...,2Nq, 



(80) 



Here ip'^^'''^\q) are arbitrary scalar functions of one scalar argument. In particular, the derived 
formulae describe the big class of solutions having the form of rational functions of exponents, 
such as solitary wave solutions. The simplest examples of them will be given below, see eqs.flSSl 



[931I96D . 

In the scalar case (corresponding to eqs.([I]) and ([2])) one has M = 1 so that eq. (J79i) reads: 

2 

i=i 

with the same expressions for /c*^*^ given by eqs. (l80|) . 

The reductions to the lower dimensional PDEs. Note that the described solution space 
admits reduction ([3]) embedded into the following set of reductions 

=^Jf = Vz,j, (82) 
where n is either 2 or 3. In fact, similar to reduction (17^ . these reductions mean 

€7''\<D - 5 ^t''""^) (83) 

in eq. (!79!) . For instance, let n = 2 which corresponds to reduction ([3]). Then we must substitute 
the following expression into eq. (!79!) : 



where ip^^'^^^ {k = 1,2) are arbitrary constants. In the scalar case (corresponding to nonlinear 
equations (jl]) and ([5])) eq.l lH^ reads 

i/j(^"'''\q) = iIj^^"'''^^ (^q - VA(«=';™)) + (^q + VaM^^ , (85) 

which must be substituted into eq. flST]) . 

The simplest solution to eqs.([T]) and ([2]). The simplest nontrivial solution to eqs.([T]) and 
([2]) corresponds to A^o = '^o = 1 and A*^"'''*") = 0, n = 2, 3. Then eq. (!8T|) reads 



2 r 
i=i ^ 



(86) 
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Eqs.( 1771) become the system of two following equations: 



Their solution reads: 



A 



^(l;ll) 



^J2;ll) ^(2;11) 



Ai 



^a;ll) ^(1;11) 
^(2;11) ^(2;ll) 



w 



(1;11) 

y 

(2 

y 



^(2;ll) ^{2;ll) 



(87) 

(88) 
(89) 



have four arbitrary functions of single variable ^jJ^'^^''^\q) , l,i = 1,2. 



Formulae 



Let us write the explicite formulae for the particular case ijj^''^'^\q) = ipii5{q — qu), where 
both ^pii and qu are arbitrary constants: 



A = ^ eQi,+Q2, _ (-l)i+iK2,)V'(^^'*V^''''\ 



2 



Ai 



A^{ii;^)^(2i;.)^ 



(90) 
(91) 

Qu = qiiX + k^'^'^^y + qlh + qlh + k^'^''\r, (92) 
fcai;i) = ^ + K^^, fcai;2) = ^ _ i^z. = ^7i^o' + 4gz. + 4AZo, ^,/ = l,2. (93) 

The solution v has no singularities if A 7^ 0, which is true if, for instance, the following relations 
are valid: 



Ki, >i^2. >0, 2 = 1,2, V^(^1^2)^Q^ ^(ll;l)^^(21;l)^^(21;2) 

For example, let 



K2i = q2i = ^, i = l,2 



4 



r, V'ii = ei>0, V'i2 = -6<0 



Then 



(94) 



(95) 



(96) 



which is the kink. 



The simplest solution to eqs.(I3j) and 1^. We take eq. (l5Tl) with A^o = '^o = 1 and 

^(3;/;m) _ g. 

v[/(';ii)(f) = J2 [ dq^^''-^'\q)e'^''^'''"'''y''^'^^^^^^^ 1 = 1,2. (97) 

i=i ^ 

Note that eqs.(jl]) and ([5]) require reduction ([3]), which means that function 7/;('^'*)(g) in eq. fl97|) 
is given by eq. flHSj) . For instance, let 

^C™;*^) = S^.iJlrm. (98) 
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Then one has expressions 

2 

A = ^e«-+«^^((-l)^+iKH-(-iy+iK2,)V'm^2ii, (99) 

instead of expressions (I90H92I) . while relations (!93l) remain the same. The conditions similar to 
(p^ with replacement tp^^^''^^ tp^jf^ guarantee that the solution v is non-singular in this case 
as well. 

Let, for example, ^212 = 0, ^/^m = ^1 > 0, -0112 = — ^2 < and Ku > K21 > 0. Then one 

has 

eQ"6(gii - g2i) + eQ^^6(gii + 921) , . 

^ eQn(i^:n-i^2i)ei + eQ-(Kn + iC2i)e2' ^ ^ 

which is the kink. 



6 Obstacles to the complete integrability of eqs.([T]) and 

m 

It is important to note that our algorithm does not describe the full solution spaces to nonlinear 
PDEs (II]), ([2]), (jlj) and ([5]). To justify this statement we consider eqs.dl]) and ([2]), while eqs.(jl]) 
and dS]) may be treated in a similar way. The simplest argument is following. By construction, 
eqs.(II]) and ([2]) must be commuting flows, i.e. the expression 

[eq.m]t, - [eq.m]yr (101) 

must be identical to zero. However, the direct calculation shows that this expression yields 
the nonlocal three-dimensional PDE having rather complicated form (we do not represent it 
here). This means that the constructed solution space to the equations and ([2]) is two- 
dimensional (in other words, it may depend on the arbitrary functions of two independent 
variables) while the full solution space to four-dimensional PDEs and ([2]) must be three 
dimensional. The dimensionality of the solution space is confirmed in Secj5]by formulae (IHTIl . 
which shows us that the solution space depends on ANqUo arbitrary functions of single variable 

m=l,...,no, z = l,2. Since both A^o and no are arbitrary positive 
integers and may go to infinity, these functions may approximate arbitrary functions of two 
variables in the solution space. Thus, the possibility to have an arbitrary function of three 
independent variables in the solution space (which would provide the full integrability) remains 
an open problem for further study. Without details, we remark that solution space to four- 
dimensional eq.dl]) (available due to our algorithm) is two-dimensional as well, while solution 
spaces to three-dimensional equations (jl]) and are one-dimensional, which may be obtained 
owing to the formulae (IHTl) and (1851) . 

Since eqs.(IT]) and ([2]) have been derived from the matrix KP and the later is well known 
(2+l)-dimensional PDE integrable by the ISTM, one can expect the same type of integrabil- 
ity of eqs.dH) and ([2]). In particular, one can expect that these equations are compatibility 
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conditions of some overdetermined linear spectral problem derivable from the overdetermined 
linear spectral problem for the matrix KP. However, this is not true. In the next subsection 
we consider eq.(l2]) as an example and show that the linear spectral problem for this equation 
is not well defined. Thus, at the moment, the ISTM is not a suitable tool for solving eq.([2l). 
The same conclusion is valid for eq.([T]), (jlj) and ([5]) as well. 

6.1 The spectral problem for the matrix KP supplemented by dif- 
ferential reduction (12111261) and Frobenious type reduction (1481) 

Although the algorithm described above allows one to find a big manifold of solutions to eqs.([T]), 
([2]), (jlj) and ([5]) it does not gives us an algorithm to derive the linear spectral problem for any 
of these equations, i.e., at the moment, we are not able to obtain such overdetermined system 
of linear equations for some spectral function whose compatibility condition results in eq.([T]), or 
eq.(I2]), or eq.(jl]), or eq.([5]) without additional requirements to the coefficients of linear system. 
As an example explaining this problem, we consider eq.([2]) as nonlinear equation derivable from 
the matrix KP f|T71) after differential reduction (introduced in Secj3]) followed by the Frobenious 
type reduction (introduced in SecH]) in this subsection. 

The obstacle to derive the linear spectral problem for eq.(l2]) is assotiated with the spectral 
representation of differential reduction ( l2T]|26i) which is a reduction of the potentials of the linear 
spectral problem for the matrix KP. In other words, we do not know which reduction must be 
imposed on the spectral function to generate the above differential reduction for the potentials 
of the spectral problem. Thus the problem of spectral representation of the eq.([2]) remains 
open as well as the problem of its complete integrability. Nevertheless, we derive some linear 
spectral problem whose compatibility condition leads to eq.(I2]) imposing the above differential 
reduction " by hand" . 

The spectral problem for matrix KP f|T7j) reads 



where ip is a spectral function and A is a spectral parameter. Differential reduction fl21|26p 
requires the following block structure of ip: 



where tp^^^'' ^ire uqM x uqM matrices. So that equations fll02p acquire the following block 
structure: 



F(2) := ^,3(A; x) - ^,,,iX; x) - 3^.(A; x)Wi'\S) + \^{X- x) {wil\x) - W^X^)) 



(102) 




(103) 




(104) 



where F^^'^"^^ are rioM x rioM matrix equations. The spectral problem is now represented by 
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the first-row blocks of eq.f ll04p : 



_p(i;ii) 
^(i;i2) 



) + + 2^(^^).« + 2^(i^« = 0, 



w 



(1) 

t2 



AD 

fxx 



(11) - „(i) 



^t2 



(12) f ^(1) _ ^(1) 



^xx 



0, 



0, 



(105) 
(106) 
(107) 

(108) 



wliile tlie second-row blocks of eq. fll04p are equivalent to the system ( ll05H108p up to the re- 
placement of the spectral functions ?/)(i^) ■?/'(^^\ A; = 1,2. Emphasize that the differential 
reduction is already imposed on the coefficients of the above linear system due to the special 
form of the functions p^^^ and g(i) given by eq. fl26|) . It may be shown by the direct calculations 
that the compatibility conditions of eqs. fll05til08l) yield the system of nonlinear PDEs which is 
equivalent to system fHS]) . fHTI) . i.e. we derive two complete compatible systems of nonlinear 
PDEs so that eqs. (HSj) describe the t2-evolution while eqs. (H7|) describe the ts-evolution of fields. 
However, since differential reduction has not spectral representation, these nonlinear systems 
may not be integrated by the ISTM. Of course, this means that the eq.(l2]) (which has been 
derived from the above system of nonlinear PDEs after the Frobenious type reduction (H5]l ) 
may not be integrated by the ISTM as well. However, for the sake of completeness, we consider 
the Frobenious type reduction fHHl) imposed on the linear system (1105111081) . 

First of all, note that this reduction requires the following structure of the functions -0(11) 
and -0(1^); 



(In) 



^(n;ll) 
^(n;nol) 



^{n;lno) 



, n = l,2, 



(109) 



nonoj 



where ■0("''*-?) are M x M functions. Eqs. (1105111080 acquire the following block structure: 



p{n;lm) 



p{nm;nono) 



n,m = 1, 2, 



(110) 



where F("'^'*-^) are M x M matrix equations. Similarly to eq. (11041) . only the first-row blocks of 
these equations represent the system of independent spectral equations: 



^(ll;ri) 
p{12;n) 

p{22;n) 



I (1;") 



+ V'i','"^ + 27/;(i;iV(i'") + 2V^^^'^^p^^'^"^ + 2V^-'-'p- 



(2;2) (l;2n) 



I (1;") 



W- 



{l;n) 



(l;n) 



wi:j.'"' I -|- -0 



(2;1) fpjlil") 



'Qx 
x fx 

',(2;2) 



0, 
0, 



{l;2n) 



„(l;2n) 
fxx 



I (2;n) 



,(2;n) 



'^xxx " \^x ^x 



3(^(^^^).(^-) + ^(^^^)g(^^-) + ^r^g(^^-)) + 



7/;( 



't2 



{2;2) / „(l;2n) _ Jl;2n) 



Tx 



(111) 

(112) 
(113) 

= 0, 

(114) 

0, 
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where n = 1, . . . , no and 



P 



= w, 



(l;fc+n2(l)) 



(l;2fc) 



,(i;2fe) 



= w 



(l;fe) 



(115) 



.(i;ifc) 



= V, 



{l;fc+l) 



= V 



{l;fe) 



Here we take into account eq. fl32l) and take ni(l) = 1 without loss of generahty. The system 
of hnear PDEs f llllllll4i) . n = 1,2 represents the overdetermined system whose compatibihty 
condition results in the complete system of nonlinear PDEs which includes equations ( l65l) . ( 1671) 
and (!68|) as a complete subsystem, whose scalar version (M = 1) results in eq.([2]). 

7 Conclusions 

We have constructed two examples of the four- dimensional nonlinear PDEs starting with the 
dressing method for the (l+l)-dimensional C-integrable Biirgers hierarchy and using combina- 
tion of the Frobenius type and the differential reductions. One of these examples, eq.([2]), has 
the same dispersion relation as the KP does and may be referred to as the KP type equation. 
As a consequence, its lower dimensional reduction, eq.([S]), has the same dispersion relation as 
the KdV does and may be referred to as the KdV type equation. Although the derived four- 
dimensional PDEs are not completely integrable by our technique (see Secj6]), we are able to 
supply a big solution manifold to them with solitary wave solutions as most simple examples. 
The new feature of this algorithm in comparison with one represented in [16] is the differential 
reduction which introduces the new independent variable y into the nonlinear PDEs. In turn, 
this variable allows one to introduce the set of new r- variables by formula fl28|) . which is a new 
method of increasing the dimensionality of solvable nonlinear PDEs. An important problem is 
to overcome the obstacle to the complete integrability of the derived four-dimensional nonlin- 
ear PDEs (|T1[2l) and their three dimensional reductions (j^fSj) . It is also important to find the 
physical application of the derived nonlinear PDEs. Regarding this problem we must note that 
the additional independent variable y appears only in the nonlinear terms of eqs.([2]) and 
which usually means the existence of solutions with break of wave profiles. So that one can 
expect that these equations describe the break of wave profiles in the physical systems where 
the KP and the KdV appear. 

The work was supported by the RFBR grants 07-01-00446, 08-01-90104 and 09-01-92439 
and by the grant NS-4887.2008.2. 
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